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Radiative decays of the ip(2S) into jKKtv and jr/TT + TT^ final states are studied using 14 million 
il>(2S) events collected with the BESII detector. Branching fractions or upper limits on the branching 
fractions of ip(2S) and \cj decays are reported. No significant signal for ?/(1405)/j/(1475) is observed 
in the KKtt or r)ir + -K~ mass spectra, and upper limits on the branching fractions of tp(2S) — > 
7??(1405)/??(1475), »j(1405)/»y(1475) -> KKtt and j/tt+tt - are determined. 

PACS numbers: 13.25.Gv, 12.38.Qk, 14.40. Gx 



I. INTRODUCTION 

ip(2S) decays via three gluons or a single direct pho- 
ton have been extensively studied [1]. However, there 
have been fewer studies of ip(2S) radiative decays [2]. 
Further study of ip{2S) radiative decays will provide 
more information about the ip(2S) decay mechanism 
and may help in understanding problems like the u pn 
puzzle". The "12% rule" predicted by perturbative 
QCD [3] is expected to be applicable to tp(2S) radia- 
tive decays [4] , so it can be tested by measuring more 
of these decays. Furthermore, if the 12% rule is obeyed 
for the ip(2S) — > 777(1440) [2] decay, we might expect 
to observe 77(1440) in tp(2S) decays into -fKKir and 

A glueball candidate, the 7/(1440), is now re- 
garded as the superposition of two independent states, 
the 77(1405) and the 77(1475), with different decay 
modes [2]. The 7/(1475) could be the first radial exci- 
tation of the 7?'(958), while the 7/(1295) could be the 
first radial excitation of the 77. The results of L3's 
measurements on the KKtt and 7/7r + 7r~ channels in 
77 collisions suggest that the r\ (1405) has a large glu- 
onic content [5]. However CLEO did not confirm L3's 
results with a five times larger data sample and set 
upper limits on r 77 (7/(1405)),8(7/(1405) -» KKtt) and 
r 77 (r/(1475))S(77(1475) -» KKtt), which arc still con- 
sistent with the glueball and the radial excitation hy- 
potheses for 77(1405) and 7/(1475) [6]. 

Many studies have been made for 7/(1405)/7/(1475) 
with J/tjj decays into KKtt, i]tt + tt~ , 4tt, and jp° [2], 
while in ip(2S) decay, only MARKI reported an 
upper limit at the 90% confidence level (C.L.) for 
V>(2S) -» 77/(1405) -» jKKtt [7]. Here we study 
t/(1405)/t/(1475) in ip(2S) radiative decays to jKKtt 
and -fr]TT + TT~ final states using a sample of 14 x 10 6 
ip(2S) events. 

In lowest-order perturbative QCD, the x c o and 
Xc2 decay via the annihilation of their constituent cc 
quarks into two gluons, followed by hadronization of 
the gluons into light mesons and baryons, so these de- 
cays are expected to be similar to those of a gg bound 
state, while xa cannot decay via the annihilation of 
their constituent cc quarks into two gluons. So sys- 
tematic and detailed studies of hadronic decays of the 



Xc.j may help in understanding the decay patterns of 
glueball states that will be helpful for their identifica- 
tion. 

BESI collaboration studied XcJ decays into 
KgK + TT~ + cc. [8] and reported Xci branching frac- 
tions and upper limits on branching fractions of XcO 
and Xci decays. In this paper, we report measure- 
ments of ip{2S) decays into jKKtt and 77/7r + 7r~ final 
states using 14 million ip(2S) events collected with the 
BESII detector. With about a four times larger ip(2S) 
sample, more precise results are expected. Branching 
fractions or upper limits of tp(2S) and XcJ decays are 
reported. 



II. THE BESII DETECTOR 

The Beijing Spectrometer (BESII) is a conventional 
cylindrical magnetic detector that is described in de- 
tail in Ref. [9]. A 12-layer vertex chamber (VC) 
surrounding the beryllium beam pipe provides input 
to the event trigger, as well as coordinate informa- 
tion. A forty-layer main drift chamber (MDC) located 
just outside the VC yields precise measurements of 
charged particle trajectories with a solid angle cov- 
erage of 85% of 47r; it also provides ionization energy 
loss (dE/dx) measurements which are used for particle 
identification. Momentum resolution of 1.7%yl + p 2 
(j> in GeV/c) and dE/dx resolution for hadron tracks 
of <~ 8% are obtained. An array of 48 scintillation 
counters surrounding the MDC measures the time of 
flight (TOF) of charged particles with a resolution of 
about 200 ps for hadrons. Outside the TOF coun- 
ters, a 12 radiation length, lead-gas barrel shower 
counter (BSC), operating in limited streamer mode, 
measures the energies of electrons and photons over 
80% of the total solid angle with an energy resolution 
of a E /E = 0.22/y/E (E in GcV). A solcnoidal mag- 
net outside the BSC provides a 0.4 T magnetic field 
in the central tracking region of the detector. Three 
double-layer muon counters instrument the magnet 
flux return and serve to identify muons with momen- 
tum greater than 500 MeV/c. They cover 68% of the 
total solid angle. 
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III. EVENT SELECTION 

The decay channels investigated in this paper are 
ip(2S) — > ■yKgK+n' + c.c, jK+K^tt" and jijit+tt^, 
where K s decays to 7r + 7r _ , 77 to 77, and 7r° to 77. 

A neutral cluster is considered to be a photon can- 
didate if the following requirements are satisfied: it 
is located within the BSC fiducial region, the energy 
deposited in the BSC is greater than 50 MeV, the first 
hit appears in the first 6 radiation lengths, the angle 
in the xy plane (perpendicular to the beam direction) 
between the cluster and the nearest charged track is 
greater than 8° , and the angle between the cluster de- 
velopment direction in the BSC and the photon emis- 
sion direction from the beam interaction point (IP) is 
less than 37°. 

Each charged track is required to be well fitted by a 
three-dimensional helix, to have a momentum trans- 
verse to the beam direction greater than 70 MeV/c, 
to originate from the IP region (V xy = \jv x + V y < 2 

cm and \ V Z \ < 20 cm) if it is not from K s decay, and to 
have a polar angle | cos#| < 0.8. Here V x , V y , and V z 
are the x, y, and z coordinates of the point of closest 
approach of the track to the beam axis. 

The TOF and dE/dx measurements for each 
charged track are used to calculate Xp/dW values 
and the corresponding confidence levels Probpip>(i) 
for the hypotheses that a track is a pion, kaon, or 
proton, where i (i — ir/K/p) is the particle type. For 
each event, charged kaon candidates are required to 
have ProbpiD (K) larger than 0.01, while charged pion 
candidates are required to have Probpio^) > 0.01. 



IV. EVENT ANALYSIS 
A. rp(2S) -> iK%K +7T" + c.c. 

For the final state "fK ± ir T TT + TT^ , the candidate 
events are required to have at least one photon candi- 
date and four good charged tracks with net charge 
zero. A four constraint (4C) kinematic fit is per- 
formed to the hypothesis ip(2S) — ► jK^irTir+n^ , and 
the x 2 01 the fit is required to be less than 15. If 
there is more than one photon, the fit is performed 
with the photon candidate which has the largest en- 
ergy deposit in the BSC. A 4C-fit to the hypoth- 
esis ip{2S) — > 77T + 7r~7r + 7r~ is also performed, and 
xlcil-K ± 7r =F 7r + 7r _ ) < x! c (77T + 7r - 7r + 7r~) is required 
to suppress background from 77r + 7r~7r + 7r~. 

Backgrounds from ip{2S) — > 7r + 7r~ J ftp are rejected 
with the requirement \m^ ml - 3. 1 1 > 0.05 GeV/c 2 , 



0.550 fTTTTTTTT-rp 



~ 0.525 



CD 

o 

X 0.500 

E 0.475 



0.450 La - 



I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 



I "I 



0.00 0.05 0.10 0.15 0.20 0.25 
Lxy (m) 



FIG. 1: The scatter plot of tx + -k invariant mass versus 
the Kg decay length. 



where w^ecDii is the mass recoiling from each possi- 
ble 7r + 7r~ pair. Figure 1 shows the scatter plot of 
7r + 7r~ invariant mass versus the decay length in the 
transverse plane (L xy ) of Kg candidates. A clear Kg 
signal is observed. Candidate events are required to 
have only one Kg candidate satisfying the require- 
ments 1^^+,,- —0.498| < 0.015 GcV/c 2 and L xy > 0.5 
cm. After Kg selection, if one of the remaining tracks 
has a momentum higher than 1.5 GeV/c, it is taken 
as a charged kaon. Otherwise, the track types are se- 
lected using their Xk„ values, i.e., if Xr+ v - < xl+ K - » 
the final state is considered to be jK s K + ir~ ; if 
Xk-^ < X^-k+i the final state is considered to be 
jK s K-7r+, where X k = Xpid(K) + XpidW- 

With this selection, Fig. 2 shows the mass dis- 
tribution of K s K + -k~ and K s K~ir + for candidate 
events. There is a clear Xci signal, but no clear 
r/(1405)/r/(1475) signal. The biggest background con- 



tamination comes from ip(2S) — > 7r° K s K + n~ + c.c, 
which is estimated with the data sample, and the 
other backgrounds are estimated by Monte Carlo 
(MC) simulation. 

In the high mass region, the fit of the K s K + n~ + 
c.c. invariant mass spectrum is performed after sub- 
tracting the known background, and a second order 
polynomial is used to describe the shape of the re- 
maining unknown background (see Fig. 3). The Xco 
peak is described with a Breit-Wigner folded with a 
double-Gaussian resolution function determined from 
MC simulation, while the Xci and Xc2 peaks are de- 
scribed only with double-Gaussians resolution func- 
tions because their widths are much smaller than the 
mass resolution. The masses of the three XcJ states 
and the width of x c0 are fixed to PDG values [2]. 

A binned maximum likelihood method is used to 
fit all events with KgK^ir 11 mass between 3.2 and 
3.65 GcV/c 2 . The numbers of events are 3.9 ± 4.6, 
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220 ± 16, and 28.4 ± 7.6 with statistical significances 
of 0.9d, 22.0cr, and 4.8(7 [10] for Xco, Xci, and Xc2, 
respectively. 
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FIG. 2: Invariant mass distributions for ip(2S) — > 
KgK + ir~ + c.c. candidate events in the low mass region 
(upper plot) and high mass region (lower). Dots with er- 
ror bars are data, and the hatched histogram is simulated 
background. 



80 



> 

CD 

o 

T^40|- 
co 

LU 



LU 




3.30 3.40 3.50 3.60 3.70 

m(KgK + jt-+c.c.) (GeV/c 2 ) 



FIG. 3: The result of the K%K + -r- + c.c. mass fit. The 
curve shows the best fit described in the text. 

Figure 4(a) shows the Dalitz plot of \d — * 
KgK + n- +c.c. candidate events with 3.48GeV/c 2 < 
m K o K±nT < 3.53GeV/c 2 . The clusters of events in- 



shows the K ± n T invariant mass distribution after an 
additional requirement m K o n ± > 1.0 GeV/c 2 to reject 
K^KT events. Fi gure 4(c) shows the K^tt invari- 
ant mass after the requirement m K ± n T > 1.0GeV/c 2 
to reject K*°K events. 

For Xci — * KgK + n~ + c.c. candidate events, 
the if ± 7r =F and K^tt^ 1 mass spectra are fitted with 
A"* (892) and AT}(1430) signal shapes determined from 
MC simulations plus a threshold function for back- 
ground. For A"*(892)° and K*(892) ± , the fitted 
numbers of events are 22.5 ± 7.3 and 26.7 ±11.0 
with corresponding statistical significances of 3.5<7 and 
3.0(7, respectively. For if}(1430), there are three 
states: K*(U30), K*(U30), and if*(1410) around 
1430 McV/c 2 . With the detection efficiencies aver- 
aged with equal weight, the numbers of events includ- 
ing these three hypotheses are calculated to be 22 ± 15 
and 45 ± 26 for if}(1430)° and if}(1430) ± , respec- 
tively [11]. The upper limits on the numbers of events 
at the 90% C.L. are calculated to be 41 and 79. 

In the low mass region, no ?7(1405)/?7(1475) sig- 
nal is observed in the KgK^n^ invariant mass dis- 
tribution. Here the fit is performed under two hy- 
potheses: one for r/(1405) with mass 1410 McV/c 2 , 
width 51 MeV/c 2 , and mass resolution 7.1 MeV/c 2 ; 
the other for ?7(1475) with mass 1476 MeV/c 2 , width 
87 MeV/c 2 , and mass resolution 7.7 MeV/c 2 . The 
KgK + Tr~ + c.c. invariant mass distribution is fitted 
with a Breit-Wigncr folded with a Gaussian resolu- 
tion and a second order polynomial for background. 
The mass, width, and mass resolution are fixed to the 
values above. The signal is very weak, so upper limits 
on the number of events are calculated to be 11 and 
16 for 77(1405) and r?(1475), respectively. 



B. ip(2S) -> -yK+R-n 

For this channel, candidate events are required to 
have two charged tracks with net charge zero and three 
photon candidates. A 4C-fit is performed under the 
ip(2S) — » ^"fK + K^ hypothesis, and the x 2 of the fit 
is required to be less than 15. The invariant mass of 
the charged kaon tracks is required to be less than 3.0 
GeV/c 2 to veto ?p(2S) — > neutral + J/tp background. 
With three selected photons, there are three possible 
combinations to reconstruct ir°, and the combination 
with invariant mass closest to m^o is taken as the 
7T° candidate. Figure 5 shows the 77 invariant mass 
distribution, where a clear n° signal is observed. 

After requiring |m 77 — m^o] < 0.03GeV/c 2 , Fig. 6 



dicate if* (892) and AT} (1430) signals. Figure 4(b) events. There is no jj(1405)/jj(1475) signal in the low 
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FIG. 4: (a) Dalitz plot of Xci — > ifs-/\ + 7r~ +cc. candidate 
events. The (b) K tt t and (c) ifgTr* invariant mass dis- 
tributions. In (b) and (c), dots with error bars are data, 
and the histograms show the best fits described in the text. 



mass region. Upper limits on the number of events 
are calculated to be 9 and 9 for 7/(1405) and 7/(1475), 
respectively. 
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FIG. 5: The 77 invariant mass distribution for tp(2S) 
^y^K + K ~ candidate events. 
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FIG. 6: The K+R-tv invariant mass distribution for 
ip{2S) — * jK + K~tt° candidate events. Dots with error 
bars are data, and the hatched histogram is the simulated 
background. 



V(2S) 



7r77T + 7T 



The final state of this channel is 7r + 7r~777. Events 
with two charged tracks with net charge zero and three 
photon candidates are selected. A 4C-fit is performed 
for the hypothesis *p(2S) — ► 7r + 7r~777, and the \ 2 
of the fit is required to be less than 15. Background 
from tp(2S) — > tt + it~ J/ip is rejected with the require- 
ment |rnrecoi/ — 3. 1 1 > 0.05 GeV/c 2 . Background from 
*p(2S) — > neutrals + J/ip is suppressed with the re- 
quirement m 77r + 7r - < 2.8 GeV/c 2 , where m 77r + 7r - is 
the invariant mass of the 7r + 7r~ and the photon which 
does not come from 77 decay. 

With the above selection, Fig. 7 shows the 77 invari- 
ant mass distribution, where 77 includes all possible 
combinations among the three photon candidates. A 
clear rj signal is observed. The smooth background 
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FIG. 7: The 77 invariant mass distribution for ip(2S) — > 
7T+7r - 777 candidate events (dots with error bars). The 
curves show the best fit described in the text. The hatched 
histogram is the 77 distribution of background events from 
the continuum and the 14M inclusive decay MC sample 
with signal events removed. 



comes from many channels and can be described by 
the sum of continuum events and ip(2S) inclusive de- 
cay MC events, where the signal events have been 
removed and some known background channels are 
replaced by MC simulated results. The main back- 
ground of the 77 signal comes from ip(2S) — > rjir + ir~ n° , 
which is estimated using tp(2S) data. We also stud- 
ied other possible channels listed in the PDG [2] that 
might contaminate the rj signal, but the contamina- 
tion is negligible. A fit of the m 77 spectrum yields 
553±60 events, and the background contamination is 
estimated to be 135±59 by fitting the hatched his- 
togram in Fig. 7. In the fit, the r\ signal is described by 
the double-Gaussian shape determined from ip(2S) — > 
jr]ir + Tr~ MC simulation. After the background con- 
tamination is subtracted, the number of r\ events is 
418±60, with a statistical significance of 7.3o\ Here 
the background contamination is subtracted from the 
total number of observed events, and the uncertainty 
on the number of background events is taken as a sys- 
tematic error. This method to deal with the back- 
ground contamination is also applied to the following 
analyses. 

An 77 candidate is defined with the requirement 
|m 77 - 0.548| < 0.04 GeV/c 2 . Figure 8 shows the 
777T+7T - invariant mass distributions in the low and 
high mass regions. Clear 77' (958) and Xci signals are 
seen. 
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FIG. 8: The rin + n~ invariant mass distributions for 
tp(2S) — > 7777r + 7r~ candidate events in the low mass re- 
gion (upper plot) and high mass region (lower). The dots 
with error bars are data, and the hatched histogram is the 
background estimated from the n sidebands. The curves 
show the best fit described in the text. 



1. ip(2S) -> 77/(958) and 7 ??(1405)/??(1475) 

Besides the r/ (958) signal, there is also a small peak 
at 1430MeV/c 2 , which could be an 7/(1405), /1 (14.20), 
p(1450), or 77(1475) listed by PDG [2]. /i(1420) and 
?7(1475) dominantly decay into KKtt, but no signifi- 
cant signal of /i(1420) or r/(1475) is observed in the 
KsKtt invariant mass distribution (see Fig. 2(a)). 
The tp(2S) — > 7p(1450) decay is forbidden by C-parity 
conservation. So the peak at 1430MeV/c 2 is assumed 
to be 77(1405) signal, and more will be discussed later. 

Assuming 7/ (958) and 7/(1405) signals, the low mass 
region is fitted with the MC distributions plus a sec- 
ond order polynomial for background (see Fig. 8). 
The fit yields 24.2±5.4 and 13.8±7.0 events, and 
the peaking background events are estimated to be 
0.9±1.4 and 4.0±4.5 from 77 sidebands. The 77 side- 
band region is defined by |m 77 — 0.38 1 < 0.04 GeV/c 2 
and 



0.72| < 0.04 GeV. After background sub- 
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FIG. 9: Dalitz plot of Xci — * ri-K + ir candidate events. 
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FIG. 10: The a^ir T invariant mass distribution for 
4>(2S) — > 7a (980) ± 7r T candidate events (dots with er- 
ror bars). The curves show the best fit described in the 
text. The hatched histogram is the m J]v + v - distribution 
of the events in ao(980) sideband region. 



traction, the numbers of y/(958) and r?(1405) events 
become 23±5 and 10±7, and the statistical signifi- 
cances are 6.6a and 1.4<r, respectively. 

Since the significance of r/(1405) is low and there is 
no clear 77(1475) signal, upper limits at the 90% C.L. 
on the numbers of events for 7,(1405) and ?7(1475) are 
calculated to be 24 and 20, respectively. 

2. V(25)^7Xd 

The fit in the high 771^+^- region yields 256±28 
Xd events (see Fig. 8), and the peaking background 
events are estimated to be 34±15 from the 77 side- 
band region. The 77 sideband region is defined by 
|m 77 - 0.38| < 0.04 GeV/c 2 and |m 77 - 0.72| < 
0.04 GeV. After the background contamination is 
subtracted, the number of x c i signal events becomes 
222±28, with 8.8c statistical significance. 

The Dalitz plot of Xci — * rj-K+ir^ candidate events 
within the Xci mass window (3.46-3.56) GeV/c 2 is 
shown in Fig. 9. The horizontal and vertical clus- 
ters with 77l r)7r around 1 GeV/c 2 correspond to Xci ~ > 
a (980)7r, and the diagonal band is Xci — * /2(1270)7r. 

The ao(980) ± 7r =F invariant mass distribution for 
events satisfying (\m m ± — 0.985| < 0.1 GeV/c 2 ) is 
shown in Fig 10. The distribution is fitted with a MC 
determined double-Gaussian function plus a second 
order polynomial for the background. The fit yields 
79±14 Xci candidate events, and the number of back- 
ground events contributing to the peak is estimated 



to be 21±11 by using a similar fit for events from the 
do (980) sideband region. The ao(980) sideband region 
is defined by \m vn ± — 1.6| < 0.3 GeV/c 2 . After sub- 
traction, the number of Xci signal events is determined 
to be 58±14, with a 4.5ct statistical significance. 

Figure 11 shows the scatter plot of to 77 versus 
m w + n - of the Xci candidate events. The /2(1270)?7 
signal region is defined by \m w + w - — 1.275| < 
0.185 GeV/c 2 and |m 77 - 0.548| < 0.04 GeV/c 2 . The 
/2(1270) sideband region is defined by \m n + n - — 0.7| < 
0.185 GeV/c 2 and \m n+n - - 1.85| < 0.185 GeV/c 2 , 
and the 77 sideband region by |m 77 — 0.38 1 < 
0.04 GeV/c 2 and |m 77 - 0.72| < 0.04 GeV (see Fig. 

11) . The central box in the figure is the signal region, 
the two boxes located above and below the signal re- 
gion are the 77 sidebands (named B^ and £? 2 ), and the 
two on the left and right of the signal region are the 
/ 2 (1270) sidebands (named Bj 2 and BjJ. The four 
boxes at the corners are used to estimate the phase 
space contribution (named B l ph , i — 1,2, 3, 4). We use 
the formula: 

\ x (Bl+Bl + B) 2 +B} 2 )-- A x (Bl h +B 2 ph +B s ph +B 4 ph ) 

to obtain the m vn + n - distribution of the events in the 
sideband regions. Figure 12 shows the ?77r + 7r~ invari- 
ant mass distributions of the signal and background 
regions, and a clear Xci signal is observed. 

The Xci signal is fitted with a double-Gaussian func- 
tion determined from MC simulation plus a second or- 
der polynomial to describe the background (see Fig. 

12) . The fit yields 65±13 events, and fitting the side- 
band region events yields 12±7 sideband background 
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FIG. 11: The scatter plot of m 77 versus m^^- for can- 
didate events. 
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FIG. 12: Invariant mass distribution of /2(1270)t7 for 
ip(2S) — ♦ 7/2(1270)77 candidate events (dots with error 
bars). The curves show the best fit described in the text. 
The hatched histogram is the m J17r + 7r - distribution of 
events in the sideband regions. 



events. After subtraction, the number of Xci signal 
events is 53±13 with a 4.8c statistical significance. 



Monte Carlo simulation is used for mass resolution 
and detection efficiency determination. In this analy- 
sis, a GEANT3 based Monte Carlo package with de- 
tailed consideration of the detector performance (such 
as dead electronic channels) is used. The consistency 
between data and Monte Carlo has been carefully 
checked in many high purity physics channels, and 
the agreement is reasonable [12]. 

For ip(2S) -» 77/(958) and ip(2S) -» 
777(1405/1475), the photons are distributed according 
to a 1 + cos 2 9 distribution. For i[)(2S) — > JXcJ 
under the assumption that the processes are pure 
El transitions [13], the photons are generated as 
1 + cos 2 9, 1 — i cos 2 9, and 1 + cos 2 9 for Xco, 
Xd, and Xc2, respectively. Multihadronic decays of 
t/(958), 77(1405/1475), and Xd are simulated using 
phase space distributions. 

In the MC simulation for ip(2S) — ► 7Xci, Xd — * 
a,Qir T , the width of a (980) is assumed to be 
75MeV/c 2 in the determination of the detection ef- 
ficiency. The uncertainty of the efficiency due to the 
uncertainty of the ao (980) width is taken as a system- 
atic error in the branching fraction of Xci — > clq1t t . 

The efficiencies for the determination of the branch- 
ing fractions of ip(2S) — > 7777T+7T - , Xd — > r)ir + ir~ and 
Xd — > KgK + n~ +c.c. are determined from a weighted 
average over the intermediate processes. 



VI. SYSTEMATIC ERRORS 

Many sources of systematic error are considered. 
Systematic errors associated with the MDC tracking, 
kinematic fitting, particle identification, and photon 
selection efficiencies are determined by comparing J/ip 
and ip(2S) data and Monte Carlo simulation for pure 
data samples, such as ip(2S) — > ir + n~J/ip. 

The uncertainties on the total number of ip(2S) 
events, the branching fractions of intermediate states, 
the ao (980) width, the detection efficiency, the back- 
ground contamination, and the fitting on the mass 
spectrum are also considered as systematic errors. Ta- 
ble I summarizes the systematic errors for all channels. 



Xci decays to a ir^ and /2(1270)?7 yield the same fi- 
nal state i]TT + ir~ . MC studies show that the sideband 
analysis described above separates the two channels 
without cross contamination. However, any interfer- 
ence effects are not taken into consideration because 
of the low statistics. 



VII. RESULTS AND DISCUSSION 

Tables II, III, and IV summarize the results for 
the channels measured in this analysis. Table II lists 
the branching fractions of t/;(2S i ) decays. To com- 
pare with the 12% rule, Table II also includes the 
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TABLE I: Summary of systematic errors (%). MDC, 4C, PID, 7eff., N^s), Int., MC, Bg., Fit. and K% rec. are for 
tracking, kinematic fit, particle identification, 7 detection efficiency, ip(2S) total number, the branching fractions of the 
intermediate states, MC statistics, background, fitting on mass spectrum and Kg reconstruction, respectively. 



Channel (tp(2S) ->) 


MDC 4C PID 


7 eff. 


N 4>(2S) 


Int. 


MC 


Bg. 


Fit. 


K° s rec. 


Total 


7XcO -*■ ^K%K + i\^ + c.c. 




8.0 


6.0 


- 


2.0 


4.0 


8.1 


1.3 


20.6 


18.0 


3.4 


30.8 


IXci — * r yK s K + n^ + c.c. 




8.0 


6.0 




2.0 


4.0 


9.5 


1.3 


0.2 


0.7 


3.4 


15.0 


7Xc2 — > ^yK^K + n~ + c.c. 




8.0 


6.0 


_ 


2.0 


4.0 


9.4 


1.3 


2.9 


3.9 


3.4 


15.7 


7X ci ^-yK*(892)°K° s 




8.0 


6.0 


_ 


2.0 


4.0 


9.5 


1.3 


- 


1.9 


3.4 


15.1 


7Xci ^7^(892)*^ 




8.0 


6.0 


_ 


2.0 


4.0 


9.5 


1.3 


- 


3.2 


3.4 


15.3 


7Xd -> 7^}(1430)°^g - yK° s K + n- 


+ c.c. 


8.0 


6.0 


_ 


2.0 


4.0 


9.5 


1.3 


- 


5.9 


3.4 


16.1 


7Xd -> 7^(1430)*^ 7^gir+7r- 


~ + c.c. 


8.0 


6.0 


_ 


2.0 


4.0 


9.5 


1.3 


- 


15.2 


3.4 


21.4 


777(1405) -> 7/^+77- + c.c. 




8.0 


6.0 


- 


2.0 


4.0 


- 


1.5 


- 


5.0 


3.4 


12.6 


7?7(1475) -> 7^A: + 7r" + c.c. 




8.0 


6.0 




2.0 


4.0 


- 


1.5 


- 


24.0 


3.4 


26.7 


7?7(1405) -> 7^+^-^° 




4.0 


4.0 


2.0 


6.0 


4.0 




1.7 




12.3 




15.6 


777(1475) -> jK+R-ti" 




4.0 


4.0 


2.0 


6.0 


4.0 




1.7 




20.0 




22.2 


7777T + 7r~ 




4.0 


4.0 


2.0 


6.0 


4.0 


0.7 


2.3 


14.0 


7.7 




18.7 


777'(958) 




4.0 


4.0 


2.0 


6.0 


4.0 


3.5 


1.6 


6.0 


2.9 




12.1 


777(1405) — * 77777+77" 




4.0 


4.0 


2.0 


6.0 


4.0 


0.7 


1.8 




10.9 




14.5 


777(1475) — > 77777+77- 




4.0 


4.0 


2.0 


6.0 


4.0 


0.7 


2.0 




10.0 




13.9 


7Xci — ► 77777+77" 




4.0 


4.0 


2.0 


6.0 


4.0 


9.5 


1.1 


6.7 


6.4 




16.3 


7Xd -> 70o(980)+7T- +c.c. 




4.0 


4.0 


2.0 


6.0 


4.0 


14.4 


1.0 


19.2 


6.2 




26.5 


7Xd -» 7/2(1270)77 




4.0 


4.0 


2.0 


6.0 


4.0 


10.0 


1.0 


13.2 


2.8 




19.3 



TABLE II: Measured branching fractions and upper limits (90% C.L.) for tp{2S) decays. Results for corresponding J/ip 
decays [2] and the ratio Q h — are a l so given. 



Channel {tp{2S) ->) 


n slg - 


e (%) ^ (2S) ^(xl0- 4 ) B 


7/ ^(xl0- 4 ) 


Bum y /o > 


+ - a 

7?7-7r 7T 


418±60 


8.69 


8.71±1.25±1.64 






+ - b 
7r77r 7r 






3.60±1.42±1.83 


39±7.3 [14] 


9.2±6.2 


777' (958) 


23±5 


7.58 


1.24±0.27±0.15 


43.1±3 


2.9±0.7 


777(1405) — » 77777 + 77" 


10±7 


5.06 


0.36±0.25±0.05 


3.0±0.5 


12±10 




< 24 


5.06 


< 1.0 


3.0±0.5 


< 33 


777(1475) — » 7777r + 7T~ 


< 20 


4.80 


< 0.83 


3.0±0.5 


< 28 


777(1405) jKKty c 


< 11 


4.54 


< 0.8 


28 ±6 


< 2.9 


777(1475) -» -fKKn c 


< 16 


4.58 


< 1.5 


28 ±6 


< 5.4 


777(1405) -> jKKn d 


< 9 


3.63 


< 1.3 


28 ±6 


< 4.6 


777(1475) yKKn d 


< 9 


3.54 


< 1.4 


28 ±6 


< 5.0 


a all processes in the ip(2S) — > 


77777+ 









6 all processes excluding 4>(2S) — > 7Xci — » 77/77+77-; 
c the decay mode is ^K%K + -k~ + c.c; 
d the decay mode is -yK + K~ n° . 



corresponding J/ip branching fractions [2], as well 
as the ratio Qh of ip(2S) to J/tp branching frac- 
tions for each channel. Decay of ip(2S) to 77777+77- 
is consistent with the 12% rule expectation within 
errors; decays of ip(2S) to 777(1405) — > 7777+77- and 
777(1475) — > 7777+77- cannot be tested because of low 
statistics; while the other modes are suppressed by a 
factor of 2 - 4. The i/>(2S) -» 777' (958) branching 
fraction with 7y'(958) — > 7777+77- is more precise than 
(2.00 ± 0.59 ± 0.29) x 10~ 4 measured by BESI [15]. 
No signal for 77(1405)/t7(1475) is observed in either 



r )K { gK + ir^ +c.c. or jK + K~ir° final states. There is a 
small peak at 1430McV/c 2 in the 77777+77- final state, 
and we have treated it as 77(1405) signal. Because of 
its low statistics, we also set the upper limit at the 
90% C.L. for ip(2S) -» 777(1405)/t/(1475) -» 77777+71--. 
As shown in Table II, the upper limits at the 90% 
C.L. on ip(2S) -» 777(1405)/7/(1475) -» ^KKtt and 
7777T + 7r~ are at the same level 0.8 <~ 2.0 x 10 -4 . 

In the above study in the high 7777+77- mass re- 
gion, only Xci is considered. If we fit m m + v - with 
Xc0,i,2 together, the fit yields —32 ± 28, 250±32, and 
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TABLE III: Branching fractions for XcJ -> K%K + -k' + c.c. and xa -> 7771+71". Here B{ip(2S) -> 7Xco) = (8.6 ± 0.7)%, 
B(ip{2S) -> 7Xci) = (8.4 ± 0.8)% and B{tp(2S) -> 7 Xc 2 ) = (6.4 ± 0.6)% are used in the calculation. 



Channel 


XcJ 




e(%) 


B(xl0" 3 ) 


BESI (xl0~ 3 ) 




XcO 


< 10 


6.24 


< 0.3 


< 0.71 


Tfg^+TT- + C.C. 


Xcl 


220 ± 16 


6.80 


4.1 ±0.3 ±0.7 


2.46 ±0.44 ±0.65 




Xc2 


28.4 ± 7.6 


5.82 


0.8 ±0.3 ±0.2 


< 1.06 




XcO 


< 32 


6.64 


< 1.2 






Xcl 


222±28 


7.90 


6.1±0.8±1.0 






Xc2 


< 48 


7.17 


< 2.2 





TABLE IV: Branching fractions of X ci -> K*K, a (980)7T and / 2 (1270)r7. Here B(ip(2S) -> 7Xci) = (8.4 ± 0.8)% is used 
in the calculation. 



Channel 


71 Sl9 ' 


e(%) 


/3(xl0- 3 ) 


K*{892)°T' + c.c. 
K*(892) + K~ +c.c. 


22.5 ±7.3 
26.7 ± 11.0 


7.67 
6.20 


1.1 ±0.4 ±0.2 
1.6 ±0.7 ±0.3 


A-}(1430)°7f ' ± c.c. -> K%K + -k~ ± c.c. 
K*,(US0) + K- + c.c. -> KIK+tt- + c.c. 


< 41 

< 79 


6.28 
5.00 


< 1.0 

< 2.5 


ao(980) + 7r~ + c.c. -> r;7r + 7r" 
h (1270)7, 


58±14 
53±13 


6.10 
6.55 


2.0±0.5±0.5 
2.1±0.5±0.4 



17±26 for XcO, Xci, and X c2, respectively. The differ- 
ence in the number of X ci events is 2.3%, which has 
been taken into account as a systematic error. The 
90% C.L. upper limits on the number of Xc o and xci 
events are calculated to be 32 and 48, and the relative 
systematic errors are 12.4% and 13.3%, respectively. 
The corresponding upper limits at the 90% C.L. on 
the branching fractions are listed in Table III. 

For the XcJ KgK + ir~ +c.c. decays (listed in Ta- 
ble III), we get higher precision results compared to 
the BESI experiment [8]. The branching fraction of 
Xci — * KgK + n~ + c.c. is consistent with the BESI 
result within la, while the results for Xc j — » rjir + Tr~ 
decays are all first measurements. The branching frac- 
tions of Xci decays into intermediate processes listed 
in Table IV are all also first observations. 

XcO is forbidden to decay into KKtt or rjir + n^ by 
spin-parity conservation, and only upper limits at the 
90% C.L. are determined for these branching frac- 
tions. For Xcj decay into hadrons in lowest-order, 



Xci decay is suppressed by a factor a s compared with 
Xc2 decay. However, the branching fractions of X ci de- 
cays into KKtt and rjiT + i:~ are both much larger than 
those of X c2 decays. This result needs explanation. 
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